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We consider the cross-correlation between free electrons and neutral hydrogen during the epoch 
of reionization. The free electrons are traced by the optical depth to reionization r while the 
neutral hydrogen can be observed through 21 cm photon emission. As expected, this correlation is 
sensitive to the detailed physics of reionization. Foremost, if reionization occurs through the merger 
of relatively large halos hosting an ionizing source, the free electrons and neutral hydrogen are anti- 
correlated for most of the reionization history. A positive contribution to the correlation can occur 
when the halos that can form an ionizing source are small. A measurement of this sign change in the 
cross-correlation would directly measure the bubble bias and as, such, the halo mass. We estimate 
the signal-to-noise of the cross-correlation using the estimator for inhomogeneous reionization rim 
proposed by Dvorkin and Smith (2009). We find that with upcoming radio interferometers and 
CMB experiments, the cross-correlation is measurable going up to multipoles i ~ 1000. We also 
derive parameter constraints and conclude that, despite the foregrounds, the cross-correlation proofs 
a complementary measurement of the EoR parameters to the 21 cm and CMB polarization auto- 
correlations expected to be observed in the coming decade. 



I. INTRODUCTION 

The epoch of reionization (EoR) is one of the least 
understood periods of cosmic history, with only limited 
observational measurements (see e.g. [1 ). The absence 
of the Gunn-Peterson trough in the spectra of quasars 
implies that reionization should have been completed 
around 2: = 6 [21 13]. The total optical depth to reion- 
ization has been measured to be r = 0.084 ±0.013 [4 . If 
reionization is assumed to be instantaneous, this would 
imply a transition redshift of z^e = H- Besides these con- 
straints, we know very little about the details of reion- 
ization, such as the typical halo mass associated with the 
first ionizing objects as well as the distribution of these 
objects inside the halos. 

The spectral mapping of neutral hydrogen in emission 
[5HTT] promises to be a new probe of the EoR. The spon- 
taneous hyperfine spin flip transition causes the emission 
of a photon with a wavelength of 21 centimeters in the 
rest frame. Applying different base filters to the observed 
emission, it is possible to map the distribution of neutral 
hydrogen in the Universe as a function of redshift. The 
auto-correlation of the observed maps is very sensitive to 
the EoR parameters. 

Additionally, cross correlating the observed maps with 
other observables could provide complementary con- 
straints on the EoR parameters. For example, the 21 
cm fluctuations are expected to be correlated with galax- 
ies [12] |13 as well as with the Cosmic Microwave Back- 
ground (CMB) fluctuations through the Doppler peak 
and the kinetic SZ effect (see e.g. [14HT8] and [19] for 
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recent simulations). Unlike the fluctuations in 21 cm, 
which are a direct representation of the underlying neu- 
tral hydrogen, no direct measurement of the electron den- 
sity at high redshifts exists. The electron density can 
be measured indirectly through its integrated effect on 
the CMB, providing us with a number (r) that tells us 
the fraction of photons affected by scattering of electrons 
along the line of sight. One can go further and recon- 
struct the inhomogeneities in the optical depth field by 
considering second order effects on the CMB due to the 
screening mechanism [20 , Thomson scattering and the 
kSZ effect (see Ref. [21 , where an estimator of the 
anisotropic optical depth field is derived and Ref. [22 
for an implementation of this estimator to WMAP 7- year 
data). In this work we consider the cross-correlation be- 
tween a reconstructed map of the inhomogeneous optical 
depth Tirn (using CMB polarization observations) and a 
neutral hydrogen map measured through the redshifted 
21 cm lines. Intuitively, these two observables are ex- 
pected to be anti-correlated on most scales. 

Besides providing complementary constraints on the 
EoR parameters, the cross-correlation between the 21 cm 
field and the CMB should in principle be less sensitive to 
the details of the foregrounds. Although current 
and upcoming experiments [25] [26] are expected to be ca- 
pable of measuring the auto-correlation of 21 cm maps, 
one very persistent nuisance in extracting the signal from 
reionization are the foregrounds. Inhomogeneities in the 
21 cm signal due to patchy reionization must be sep- 
arated from fluctuations in foreground sources. Typi- 
cal foreground sources are faint radio galaxies, starburst 
galaxies and galaxies responsible for reionization. In ad- 
dition, our own Galaxy is very bright at the frequencies 
one aims at for mapping the 21 cm signal from reioniza- 
tion, exceeding the 21 cm reionization signal by several 
orders of magnitude. Attempts have been made to char- 
acterize these foregrounds [27ll2H]- Despite these efforts. 
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foregrounds can never be fully removed, simply because 
we do not know their exact origin. 

This paper is organized as follows. We review the 
physics of reionization and derive the expressions for the 
fluctuations in 21 cm aj^^ and fluctuations in the optical 
depth T£m in The former is proportional to the neu- 
tral hydrogen fraction, while the latter is proportional 
to the free electron fraction. Using a simple reionization 
model [lOl |29] where regions of HII are represented by 
spherical bubbles of typical size R we derive an expres- 



sion for the cross-correlation {rimdj^) in section In 



^IV and qV^we study the 1-bubble and 2-bubble contri- 



butions to the power spectrum. We compute the angular 
power spectrum of the cross-correlation in ^VT) We assess 
the observability of the cross correlation by using a red- 
shift weighting to maximize the signal to noise. We end 
this section with an estimate of the EoR parameter con- 
straints, when considering LOFAR and SKA noise levels. 



We present our conclusions in ^VIII In the Appendix 
we discuss the dependence of the cross-correlation on the 
parameters of the reionization model. 

Unless speciflcally mentioned, we use the following 
set of parameter values throughout the paper: h = 
0.704,116 = 0.044,11c = 0.23,11k = O.Us = 0.96 and 
r = 0.084. Lensing is included and we use the non-linear 
halo flt model to determine the power spectrum of density 
fluctuations. We use the WMAP pivot scale h = 0.002 
Mpc-i and As = 2.46 x 10"^ [4 . 



II. 21 CM BRIGHTNESS AND THE OPTICAL 
DEPTH TO REIONIZATION 

In this section we will review the standard results for 
fluctuations in the 21 cm temperature brightness and re- 
late those to fluctuations in the neutral hydrogen frac- 
tion. In the second half of this section, we will derive 
the fluctuations in the optical depth r caused by fluctu- 
ations in the free electron fraction along the line of sight, 
conflrming the results flrst obtained in [18 . 

The optical depth of a region of the IGM in the hyper- 
flne transition is given by [30 
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where = 1420.4 MHz is the rest frame hyperflne 21 cm 
{uq = A2i/c) transition frequency, Aio = 2.85 x 10~"'^^s~"'^ 
is the spontaneous emission coefficient for this transition, 
Ts is the spin temperature of the IGM, weighting the 
relative population of the atoms in the singlet state to 
atoms in the triplet state [30 , uhi is the neutral hydro- 
gen density, and v\\ the proper velocity along the line of 
sight. At high redshifts, where peculiar motions along 
the line of sight are small compared to the Hubble ffow, 
dv\\/dr\\ = H{z)/{l-\-z). At 2: = 10 dark energy and radi- 
ation are both unimportant and we can solve for H{z) in 
a matter dominated Universe, H{z) i^o^m^(l + z)^^'^. 



We can now write the following expression for the op- 
tical depth: 
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2.73(1 + z)K, = {pb-Pb)/Pb and 
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where xh is the neutral hydrogen fraction, i.e. xh = 
^Hi /{^Hi -^^e)^ pm is the matter energy density and nip 
is the proton mass. The factor (1 — 1^) addresses the fact 
that not all protons are in hydrogen, but a fraction is in 
Helium. 

The intensity along the line of sight from a thermal 
source is given by 



Jt 



(3) 



with Kj, the absorption coefficient and rjj, the emissiv- 
ity of photons. Using dl = rjydl — Hyldl = in the 
Ray leigh- Jeans limit, we have / = 2kTbiy^/c^, while 
Vu/i^u = 2kTsy'^ /c^ and /q = 2kTcmb^'^ / • Hence, we 
can write the 21 cm brightness temperature as: 

T5 = T,^be-^'' + Tsd - e-^^^) (4) 

The brightness temperature increment is deflned at an 
observed frequency v corresponding to a redshift 1 + z = 
vq/v as 
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Using Eq. we can re- write Eq. ([5| as [5l[7] 
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There are usually two types of fllters associated with 
the resolution of the experiment. First, there is a flnite 
angular resolution, which will affect all modes perpen- 
dicular to the line of sight. Second, since the brightness 
temperature of the 21 cm emission is a 3-dimensional 
fleld, we are conflned to a frequency resolution or band- 
width, which affects the modes along the line of sight (or, 
equivalently, a redshift resolution). 

The total integrated 21 cm surface brightness is given 

by 

Tbin,x) = Mx) j rfx'W^x(x')^(n,x') (7) 
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Here is an experimental band filter that is due to 
the finite frequency resolution of the instrument, which 
is centered around x (comoving distance). We define the 
dimensionless brightness temperature ip as 



Ts-T, 



cmb 



In the limit of Tg ^ Tcmhi = (1 + h)xH- 
Now To{z) can be written as 
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We will now consider fluctuations in the free electron 
density, which in turn will induce fluctuations in the op- 
tical depth. The optical depth to distance x along the 
line of sight is given by 

r{n,x) = cft I dxne{h,x)ci{x), (10) 
Jo 

where (Jt is the Thomson cross-section, Ue is the electron 
number density and a is the scale factor. Relating the free 
electron density to the free electron fraction Xg, we can 
write 



ne{n,x) 
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assuming that Helium is singly ionized. 

The average baryon density diffuses in an expanding 
background as a~^, and the free electron density becomes 



4 nip 

The optical depth can in turn be written as 
3_ .Pb,o 
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X (l + ^5(n,x0) (13) 

Therefore, we can relate fluctuations in the optical 
depth dr to fluctuations in the 21 cm brightness tem- 
perature STij [18 as 
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III. CORRELATING X AND 

We will now cross-correlate the optical depth fluctu- 
ations with the temperature brightness. As we saw be- 
fore, the CMB optical depth is proportional to the free 
electron density. If reionization is inhomogeneous, the 
free electron density is a function of position in the sky. 
Anisotropics in the optical depth produce three effects in 
the CMB: (i) screening of the temperature and polariza- 
tion fluctuations that we observe today by an overall fac- 
tor of e~^^'^\ This effect generates CMB B-mode polar- 
ization; (ii) Thomson scattering: new polarization is gen- 
erated by scattering of the local temperature quadrupole 
that each electron sees along the line of sight. This ef- 
fect also produces B-modes; and (iii) new temperature 
anisotropy is generated from the radial motion of ion- 
ized bubbles relative to the observer (the kinetic Sunyaev 
Zel'dovich effect). 

The two-point correlation function between the E- 
modes and the B-modes generated from patchy reioniza- 
tion is proportional to the anisotropic part of the optical 
depth. This fact allowed the authors in Ref. [21 to write 
a minimum variance quadratic estimator r£rn for the field 
r(n). In this work, we will use the CMB polarization fluc- 
tuations to reconstruct a map of r and cross-correlate it 
with the 21 cm field. 

We will use the shorthand notation X(n,x) = Xe{l + 
Sb). We can write the cross-correlation between the field 
X (measured through the CMB) and the field i/j (mea- 
sured through 21 cm) as 

^{xh - xjj)^S6 (15) 

Here we defined = {xh{xi)xh{x2)) - xjj, ^S5 = 
{Sb{xi)Sb{x2)) and = {h{xi)xH{x2)) . We make 
the simplistic assumption that the connected part of 
{dx^x^h^h) vanishes. Here 5x corresponds to fluctua- 
tions in the neutral hydrogen fraction, which is given by 
xh = xh{1^ 5x). 

Before we can compute ^x^^ we need to specify our 
model of reionization. We will assume that the Universe 
reionized through the growth of ionized bubbles associ- 
ated with massive halos. The bubbles themselves contain 
a single source and we assume their size to be larger than 
the non-linear scale. 

We will adopt the following average reionization frac- 
tion as a function of redshift 



where we have assumed a delta window function. 

It is worth noticing that the above expression for the 
fluctuations in the 21 cm brightness temperature is only 
valid for Tg > Tcmb- Early on, when the number of ion- 
izing sources are rare and the temperature of the IGM 
close to these sources is coupled to the kinetic tempera- 
ture by Lya photons associated with these local sources, 
this assumption breaks down, and the 21 cm signal can 
appear in absorption. We neglect this effect in this paper. 
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which is the one used in the code CAMB ^31j. Here 
y{z) — (1 + z)^/^, ^re = ^(^re) and A?/ are free parameters 
that satisfy our integrated optical depth along the line of 
sight r = 0.084. 

The ionized bubble around a given source is assumed 
to be spherical with an average radius R. We will assume 



4 



that the typical ionized bubble radii are log-normal dis- 
tributed [32 , i.e. there is a skewness towards smaller 
bubble sizes, 



where ain r is the variance of the distribution. 
The average bubble volume is then given by 

(H) = J dRP{R)Vk{R) = ^^e^<nl^ (18) 

Hence, we can define a volume weighted radius, such 
that = 47ri?o/35 which can be written as 



(19) 



If we assume that a given point in space is ionized with 
Poisson probability, we can write the ionization fraction 
as 



(20) 



with n\) the number density of bubbles. The brackets 
around ^^7g are placed to remind us that we are consid- 
ering a Poisson distribution of sources, and the result is 
averaged over the Poisson process. We further assume 
that the number density of bubbles traces the large-scale 
structure with some bias h\ 
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FIG. 1: The average ionization fraction Xe(^) as a function of 
redshift z. 



and 



dRP{R) [Vb{R)WR{kR)]\26) 



Throughout this paper we will consider bias 6 = 6 and 
Ay = 6, unless otherwise specifically mentioned. For a 
qualitative analysis of the correlation function as a func- 
tion of the parameters of our reionization model, we refer 
the reader to the Appendix. The average ionization frac- 
tion as a function of redshift is shown in Fig. [l} where 
the choice of parameters corresponds to a scenario with 
a neutral Universe at z > 13 and a completely ionized 
Universe at z < 6. 



nb{x) = n5 (1 + bSw{x)) , 



(21) 



while the average bubble number density is related to the 
mean ionization fraction as 



rib 
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The bias can in principle also be a function of R (see e.g. 
Ref. [29 ), but for simplicity we will use a fixed value 
throughout this paper. We will investigate the conse- 
quence of changing the bubble bias where it is relevant. 
6w is the matter over-density S smoothed by a top hat 
window of radius i?. 



In momentum space. 
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[sm{kR) - kR cos{kR)] , (24) 



which is the Fourier transform of Wr{x) = V^^ for x < R 
and Wr{x) = otherwise. 
We further define 



i){k) = r dRP{R)Vb{R)WR{kR) (25) 
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IV. TWO BUBBLE CORRELATIONS 

The cross-correlation between the neutral hydrogen 
and free electron fraction has two main contributions. 
First, the correlation is set by the Poisson distribution of 
ionizing sources inside the bubbles. This term is referred 
to as the 1-bubble term, and it is dominated by the shot 
noise. Second, cross correlations can also be induced by 
the enhanced probability of bubble formation (or ionizing 
sources) inside overdense regions (with probability (xg)). 
This term is referred to as the 2-bubble term, and is rel- 
evant for scales much larger than the average size of a 
bubble. 

In order to compute the 2-bubble correlation function, 
we relate the number density of bubbles to the overden- 
sities in dark matter. We define the bubble bias, which 
is related to the familiar halo bias via (see e.g. Ref. 

m) 



1 r 

fib J A 



dM ^ drih 



where drih/dlnM represents the halo mass function. 

The integral runs over all masses with some threshold 
mass scale A, above which halos are capable to form ion- 
izing sources producing a bubble with radius R. Using 
Eqs. (20) and (21) we can Taylor expand (xe) around 
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small over densities, to find 



1-(1-Xe)[l 



ln(l - xe)] + 0{6^) 



(28) 



Taking the Fourier transform of ^xx and ^^5^ we obtain 

P^',{k) = [xHHxH)b{WR){kfPss{k) (29) 
= XHHxH)b{WR){k)Pss{k), (30) 

Note that here we are considering correlations of the neu- 
tral hydrogen fraction (perturbing the free electron frac- 
tion accounts for a minus sign in Eq. (30)). The super- 
script "26" denotes the 2-bubble contribution. Also, we 
assume that the baryon fluctuations (the gas) trace the 
dark matter fluctuations. The total 2-bubble contribu- 
tion to the power spectrum of Xip results in 



Px4k) ^ -x%[\nxHb{WR){k)^l]^Pss{k) 
XH[\nxHb{WR){k) ^l]Pss{k) 



(31) 
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FIG. 2: The effective bias 6eff = Q(l - Q) at z = 12 
as a function of kR, where Q = XH^T^XHb{WR) + xh and 
Pxii^ik) = befiPss{k). The sign of the large scale correlation 
function is set by the bias (and the value of xh)- A measure- 
ment of a positive correlation at high redshift would corre- 
spond to a small bubble bias. As reionization proceeds, and 
Xh 0, Px^^jj becomes anti-correlated on large scales. 



As foreseen, the 2-bubble term is proportional to the 
total matter power spectrum with a scale dependent ef- 
fective bias 6eff- When the universe is neutral, xh 
1, the cross-correlation vanishes, Pj^^ik) = {—x'jj + 
XH)P55{k) 0, while for xh all bubbles have 
merged and P^^ 0. 

Let us further investigate the behavior of the 2-bubble 
contribution. We can write Px^{k) = Q{1 — Q)P66{k) = 
bei^Pddik)^ with Q{k^ z) = xh hiXEbiyVR) + xh- We can 
distinguish two limiting cases: for Q > 1 the power spec- 
trum effective bias 6eff is negative, representing an anti- 
correlation, while for Q < 1, ^eff is positive, and the 
2-bubble term is positively correlated. In Fig. [2] we show 
6eff for various values of the bubble bias 6 at z = 12. As 
reionization begins and ends, the 2-bubble term follows 
the matter power spectrum up to a scale-independent fac- 
tor. The transition is around kR 1, i.e. when the co- 
moving wavelength is similar to the average bubble size. 

In the large scale limit, when kR <C 1, 6eff ^ 
xniblnxH + l)(^e — blnxn)' This function changes sign 
when Xh > . In other words, the correlation func- 
tion changes sign at large scales and at early times when 
the bias b also changes, as expected. In a merging sce- 
nario, the bubble radii are small at early times and this 
should give a smaller bias. Early on, when xu ~ 1, the 
large scales are always positively correlated, because even 
if the Universe contains relatively large bubbles, the fluc- 
tuations in the ionization fraction are small compared to 
the effectively homogeneous neutral fraction, and the 2- 
bubble term traces the matter power spectrum up to a 
scale-indepedent factor. 

At scales that are smaller than the radius of the bub- 
bles, the 2-bubble term is ill-defined: the correlation 
length becomes shorter than the size of the bubbles, ef- 
fectively rendering them to one bubble. Therefore, we 
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FIG. 3: The 2-bubble term or = 1 Mpc, a bubble bias 
6 = 6, and different values of the width of the log- normal 
distribution. The relative contribution of the 2-bubble term 
to the total correlation increases rapidly with decreasing R 
and cTinR, as derived in the Appendix. 



will neglect the 2-bubble correlation term at those scales. 
In practice, we apply a smoothing filter that effectively 
cuts the correlation for /ci?o < 3. We show the 2-bubble 
contribution to the cross-correlation in Fig. [3j 



V. ONE BUBBLE CORRELATIONS 

For scales much smaller than the average bubble ra- 
dius, the correlation is dominated by the presence (or 
absence) of a single bubble [29j . The correlation between 
two (ionized) points separated by xyi = \x\ — X2\ can be 
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written as |8l [9] 

{Xe{xi)Xe{x2)) = xl ^ {Xe - xl) f {xi2 / R) , (32) 

where f{x) is a function with the following limits: 
f{x) ^ 1 for X < 1 and f{x) ^ for x > 1. If the 
probability for finding one point inside an ionized bubble 
is Xg, then when Xi2 <C R the probability of finding the 
second point in the same bubble is 1, hence their joint 
correlation probability is Xg. For large separations, the 
probability of finding two points in separate bubbles is 



the product of both probabilities, i.e., x^. Eq. (32) effec- 
tively encodes the smooth transition between these two 
regimes. The 1-bubble correlation for free electrons (or 
equivalently neutral hydrogen) then becomes: 



lb 



{Xe{xi)Xe{x2)) - xl = (Xg - xl)f{xi2/R) (33) 



As long as the bubbles do not overlap, the function 
/ can be described by the convolution of two top hat 
window functions (W^). The 1-bubble contribution to 
the power spectrum can then be written as 



pl6 



= -(Xe-Xl) miWlm + Pssik) , (34) 



where 



Pssik) = (H) 



{Wl){k')Pss{\k-k'\) (35) 



The first term in Eq. (34) is the shot noise of the bub- 



bles, which is a direct consequence of randomly placing 
ionizing galaxies in the Universe. We will later see that 
this term typically dominates the total correlation func- 
tion at late times. This can be understood by realizing 
that the bubbles tend to be larger at late times, assuming 
that the bubble size increases over time through bubble 
merging. 

Note that when correlating the free electrons with the 
neutral hydrogen, the 1-bubble contribution is always 
negative, i.e. these are fully anti-correlated when con- 
sidering just single bubbles in the Universe. In the previ- 
ous section we have also shown that the 2-bubble cross- 
correlation can be positive on the largest scales if the 
bubble bias is small. As reionization proceeds and the 
neutral hydrogen fraction decreases, the correlation on 
all scales will become anti-correlated. More importantly, 
the signal is proportional to the matter power spectrum, 
which grows as (1 + 2;)^ during matter domination. The 
bubble bias is expected to grow as a function of the bub- 
ble radius and, therefore, we expect that any positive con- 
tribution to the correlation function would only appear 
at the earliest times, when the amplitude is suppressed. 

It was shown by Ref. [29 that Pss{k) can be approxi- 
mated as 



Pssik) 



Pss{k){V,){a%) 



((H>K))2]i/2' 



(36) 



^ ^ mrj — I I 1 1 1 iii| — I I 1 1 1 iii| — I I 1 1 1 iii| — J.. I 1 1 1 iii| — I I 1 1 1 iii| 

0.01 r 

0.001 r 

"t: 0.0001 t- /' 

C\2 



cu 

%i io-« |- 

10-^ r 
10-8 
10- 



- <^in R = 0-2 ^ 

R = 0-5 = 

(7,„, = 1.0 ^ 

z = 11; R = 1 Mpc 
Liiil I I 1 1 mil I I 1 1 mil I II I mil i i 1 1 mil "1 



0.001 0.01 0.1 1 
k [Mpc-i] 



10 100 



FIG. 4: The shot noise Pf^ = {xe - xl){Vt){W^) for R = 
1 Mpc and different values of the width of the log-normal 
distribution. 
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FIG. 5: Pss as defined in the text for ^ = 1 Mpc and different 
values of the width of the log-normal distribution. This figure 
shows that this term is only relevant at small scales and does 
not contribute to the peak of the total correlation function. 



which is derived by equating the small and large scale 
limits of Eq. (35) with 



(4) 
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k^dk 



{Wl){k)Pss{k) 



(37) 



We have found the simple fitting solution of Eq. ( |36| ) to be 
accurate to the percent level for most values of {R^ crini?}- 
In the Appendix we will show that the contribution from 
Pssik) to the 1-bubble peak is relatively small for all 
parameter values in the range of interest for the r — 21 
cm cross-correlation, but is non-negligible and relevant 
at small scales. 

We show the shot noise and Pss in Figs, [i] and [sj The 
sum of these terms gives the 1-bubble power spectrum 
shown in Fig. |6] Note that the 1-bubble term from all 
three possible correlations (XX, t/^t/^ and X^Ij) is equiva- 
lent up to a sign. 

In the small scale limit, Pss{k) = Pssik) ^ and the 1- 
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FIG. 6: The 1-bubble term for = 1 Mpc and different 
values of the width of the log-normal distribution. The 1- 
bubble term has the same shape for correlations XX, Xip 
and ipip. 
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FIG. 7: k^Px^p{k)/27T^ for = 1 Mpc, 6 = 6 and ainR = 0.5. 



bubble term becomes: Px^{k) = —{xH — x'jj)Pss{k), ren- 
dering the total cross-correlation negative at these scales 
(at these scales we are applying a smoothing filter to the 
2-bubble term, so it effectively does not contribute to 
the total cross-correlation). In Fig. [t] we show the to- 
tal cross-correlation for ^ = 1 Mpc and ai^R = 0.5 at 
z = 11 (changing the redshift of the cross-correlation will 
predominantly affect the average ionization fraction Xg). 



VI. PROJECTED CROSS-CORRELATION 

Fourier transforming the dimensionless brightness tem- 
perature tp, we can write the spherical harmonic coeffi- 
cient for the 21 cm fluctuation as 



aji = M-iYj 



—^{k)af{k,z)Y;^{k), (38) 



af{k,z) = To{z) / dx'W^^,){xlMkx') (39) 
Jo 

Note that the response function is centered around 
x{z) = x^ arid in practice we take this distance to be 
somewhere between z = and z = 30. 

We can do the same for the optical depth to reioniza- 
tion, i.e. : 



—X{k)aJ{k)Y,Uk), 



with 



a^k) = (1 - Y,)aT^ T* %Mkx'), 



(40) 



(41) 



where x* corresponds to the distance to last scattering. 
Cross-correlating the two maps yields 

dk 

— ^\^{k)a}{k)af{k,z) 



(42) 



Here, 



k^Px^/i27r^). 
Let us consider the cross-correlation in the Limber ap- 
proximation. Under this approximation, we can assume 
that the Bessel functions are small, j£{x) <C 1, for x < ^ 
and peak when x ^ i. The integral over comoving mo- 
mentum k will get most of its contribution from modes 
k ^ i/X' Therefore we can make the approximation that 



^x^i^) ^ ^x^i^/x) and re- write Eq. (|42)) as 



Cj^'\z) = (i-r,) 



To{z)pb,ocrT 



dz[ 



dz[ 



Hiz[)' Jo H{z[) 

/dk 
-A%{k)je{kx{z'J)je{kx{z'2)) 

(43) 

Again, in practice we take the window function to be 
centered around < z < 30. 

In the Limber approximation we can perform the k 
integral over the product of Bessel functions as: 



f 

Jo 



dkk'^ji{kx{zi))ji{kx{z2)) 



TT 6{x{zi) - x{z2)) 



{44) 



Thus, we can write the angular cross spectrum as 

^W^ixiz')) 



c\^^\z) = (i-r,)^^^")^^'^^^ 
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FIG. 8: The cross-correlation between 21 cm temperature 
brightness fluctuations and the optical depth r at 2; = 11 for 
a model with ^ = 1 Mpc and 6 = 6. Here we used a Gaussian 
window function with bandwidth frequency Aiy = 0.2 MHz. 



We will assume that the window function is a Gaussian 
centered around redshift z with width Sx given by 



O.lMHz ) V 10 



1/2 



0.15 



-1/2 



Mpc (46) 



where Au is the bandwidth frequency of the instrument. 
We have taken into account that the power spectrum 
explicitly depends on redshift. 



We show the angular cross-correlation for Au 
MHz and several values of ain r in Fig. |8] 



0.2 



In the previous sections we have shown that the cross- 
correlation between free electrons and neutral hydrogen 
has a strong dependence on the parameters that deter- 
mine the bubble distribution as well as its bias. The 
location of the peak is set by an effective scale, which 
we derive in the Appendix. We will see that for a log- 
normal bubble distribution, this effective scale is expo- 
nential in the width of the distribution and inversely 
proportional to the average bubble radius. Therefore, 
a small change in the width of the distribution can pro- 
duce a large change in the location of the peak. We 
have also shown that at early times the 2-bubble term 
can be positively correlated at large scales. The positive 
contribution to the correlation function at large scales 
eventually vanishes when the universe further reionizes. 
However, if the bubbles are small enough, a positive con- 
tribution to the correlation function could persist until 
late times. Oppositely, if bubbles are relatively large (a 
few Mpc), the shot noise, which is negative for all scales, 
will dominate the correlation function. 



VII. IS THE CROSS-CORRELATION 
DETECTABLE? 

A. Signal-to-noise 

In this section we will determine if the cross-correlation 
is detectable. An important issue that we will address 
here are the foregrounds. As previously mentioned, the 
21 cm emission should be swamped by foregrounds, dom- 
inated by polarized Galactic synchrotron, a signal that 
is at least 4 orders of magnitude larger than the 21 cm 
brightness. We do not know the spectral dependence of 
all these foregrounds, but in general we can assume that 
they are relatively smooth in frequency along the line 
of sight, since they typically come from the same source 
(e.g. our own Galaxy). In principle, one can therefore 
remove a large part of the foregrounds by removing the 
largest modes along the line of sight (see e.g. [33 for a 
recent discussion). 

However, when cross-correlating the 21 cm field with 
the optical depth, we want to keep the largest modes, 
to which the integrated optical depth is most sensitive. 
Hence, we will keep the foregrounds in the observed maps 
and show that the cross-correlation between foregrounds 
in the 21 cm field and in the CMB should be small. In 
order to neglect the cross-correlation of the foregrounds 
between r and 21 cm, we typically need the foreground of 
the CMB to be < 10"^ times the signal [27 (given that 
^ini^ ~ lO^a^^). The synchrotron emission is roughly 
equal to the CMB signal at 1 GHz. Therefore, if we 
assume that the synchrotron scales as [34 , at 94 GHz 



synchrotron 
Im 



10"^a£5^. Thus, we 



(yV band) we have d 
estimate that the signal will be larger than the remaining 
foregrounds after cross-correlating the two maps. 

Additionally, by not removing the foregrounds, the 
21cm foregrounds will effectively act as noise term in the 
cross-correlation. In other words, even in the absence 
of correlation between foregrounds, there is still a finite 
probability that any given data point in the r map will 
correlate with a foreground measurement from 21 cm, i.e. 
the induced noise contains a term (T^m^^J/'^/), where the 
latter is the spherical harmonic coefficient of the 21 cm 
foreground map. 

Unfortunately, we do not know exactly what the level 
of synchrotron foreground is, but typically ^ M~^^ 
with 2 < a < 4. We will assume that the synchrotron 
emission scales as In Ref. [35] it was shown that 

synchrotron emission at 480 MHz has a normalized am- 
plitude of 100 mK^ < < 10000 ^^^^ with the 
actual amplitude and slope depending on the position in 
the sky. 

Although we cannot remove the foregrounds through 
implementing a large scale cutoff, we can alternatively 
try to remove a substantial part of galactic foreground 
emission. If there is a large correlation between different 
frequencies of the foreground maps, one could measure 
the foregrounds at high frequency (corresponding to a 
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completely ionized universe and, hence, with no signal in 
the cross-correlation), extrapolate with an appropriate 
scaling ^ and subtract those from the high redshift 
maps. If the correlation between different maps at high 
frequencies is of order 0.9 — 0.99, one could reduce the 
overall amplitude of the foreground by a factor of 10 — 100 
and the power by a factor of a 100 — 10^ ^. In addition, 
in Ref. [28 it was shown that down weighting the most 
heavily contaminated regions in the sky can reduce the 
effective foreground as much as a factor of 2. Note that 
this approach is different from the usual spectral fitting 
techniques l36H38j . 

We will consider a case in which the angular power 
spectrum of the foreground is given by: 



10^ 



1 



480 MHz 



(47) 



where f{z) corresponds to the frequency of the redshift 
considered and 100 < Cf < 10^ is the foreground re- 
duction factor that we can hope to achieve through a 
measurement at low redshift. 

We will assume a noise power spectrum given by [6l 
[I5l|39l[40], 



21,21 



27r 



(20 mK)2 



MHz 100 hr 

Au tint 
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Ae 




10 



9.2 



(48) 



We will do forecasts for a total integration time of 1000 
hours, and a beam with an angular diameter of A0 = 9 
arcmin. We set the bandwidth to Au = 0.2 MHz. For a 
LOFAR type experiment, we use Aeff = 10^ m^ and for a 
Square Kilometer Array (SKA) type experiment we use 
Aeff = 10^ m^. 

On the CMB side, experiments are rapidly improving 
[4TH44] , with high sensitivity experiments coming soon 
(Planck, ACTPol, SPTPol, CMBPol) and we should have 
observations of the E- and B-modes polarization spec- 
tra up to small scales in the near future. We will now 
consider a next generation polarization experiment that 
allows us to reconstruct a map of the optical depth r^m 
with the estimator proposed by Ref. [21 . This estima- 
tor was built to extract the inhomogeneous reionization 
signal from future high-sensitivity measurements of the 
cosmic microwave background temperature and polariza- 
tion fields. The authors in Ref. [21 wrote a minimum 
variance quadratic estimator for the modes of the optical 
depth field given by: 



^ Note that this is a very crude estimate. For example, it might 
be relevant to consider 21 cm signals after reionization (at low 
z) due to residual neutral hydrogen, primarily in Damped Lya 
absorbers (DLA's). 



h i2 i 

mi 1712 ^ 



X 



(49) 



where Cf^ and Cf^ are the E- and B-mode polarization 
power spectra. Nf^ and Nf^ correspond to the CMB 
noise power spectra, and are given by: 



EE 



BB 



Ap exp 



i{i + 1)^fwhm A 
81n(2) J 



(50) 



where Ap is the detector noise and 6>fwhm is the beam 
size. 

^EB 



The coupling F^^^^^ can be written as 

^EqEi 



^EB 



2i 



(2£i + l)(2£2 + l)(2^ + l) 



-2 2 



47r 
2-2 



(51) 



Here Cf^^^ is the cross-power spectrum between the 
CMB EQ-m.ode polarization without patchy reionization 
and the response field to r fluctuations Ei. [Cf^'^'^ is 
positive at large scales due to Thomson scattering, and 
negative at small scales due to the screening). 

Furthermore, the reconstruction noise power spectrum 
is given by: 



NJ 



^EB 



21- 



(52) 



We note that the main source of contamination in 
reconstructing T{n) comes from the non-Gaussian sig- 
nal from gravitational lensing of the CMB. In principle, 
unbiased estimators that simultaneously reconstruct the 
inhomogeneous reionization signal and the gravitational 
potential can be constructed [45]. For purposes of sim- 
plicity, we will estimate our cross-correlation using the 
estimator given by Eq. (49), but the results in this work 
are straightforward to generalize. 

Given that the r map is not sensitive to redshift (Ref. 
[2T] showed that the estimator is only sensitive to one 
principal component in redshift), while the 21 cm map 
can be reconstructed on redshift slices, we will give a 
weight to the cross-correlation. This weight will be built 
in order to maximize the signal to noise, in the same 
spirit as in Ref. [46 , where a weight was derived for the 
cross-correlations between CMB and Galaxy surveys. We 
write the weighted 21 cm maps as 
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(53) 



10 



10^ 
1000 
100 
10 

1 

0.1 
0.01 
0.001 
0.0001 
10-5 
10-6 
10-"^ 
10-8 
10-9 
10-10 

10-11 



I I I I 1 1 1 I I I I I I 1 1 



I I I 1 1 1 1 I I I I I I I 



10 



100 

I 




j1_ 



1000 



104 




FIG. 9: The different power spectra used to compute the 
signal to noise in Fig. |1Q| The sign oiCj^~^ has been inverted 
for the sake of comparison. The spectra are shown at z = 11. 



We then want to maximize 



and, in doing so, we find: 



wt{z) 



(54) 



(55) 



(Cf.21+iv21.21+C'/)(z)' 

which is nothing else then the projected signal over the 
projected noise. Note that we have included the fore- 
ground as a source of noise as explained before. 

We can now compute the signal to noise for the r-21 
cm cross-correlation as 

2 



TV 



/sky ^(2^+1) 



dz 



(CJ'" + iVf-)(Cf -'^ + iVf + Cl){z) 

(56) 



In Fig. 



10] we assess the level of detectability for a 

0.5 (see Fig. [o] for 



model with R = 5 Mpc and ai^R 

the angular spectra of the model considered). Again, to 
generate the spectra we use a Gaussian window function 
with Au = 0.2 MHz. We find that the signal to noise 
reaches f'^J^S/N = 8 at Cax = 3000 for a LOFAR type 
experiment, while for a SKA type experiment it reaches 

]^ /2 

fsky 'S'/A/' = 26, with /sky being the fraction of the sky 
covered. We note that for ^max ^ 1000, the signal to noise 
does not vary substantially when considering different 
values of the parameter c/, that represents the level of 
foreground subtraction. 



B. Reionization parameters 

In this section we will assess what we can learn about 
reionization by studying the 21 cm-r cross-correlation. 



FIG. 10: Total signal-to- noise for the 21 cm-r cross- 
correlation as a function of imax for a model with R — 5 
Mpc and ainR = 0.5. We consider an experiment with CMB 
noise power spectra given by Ap = 0.3/iK-arcmin and beam 
size Ofwhm = 1 arcmin. The foreground angular power 
spectrum is given by Eq. (47). We show forecasts for a SKA 
type experiment (in red lines) and a LOFAR type experiment 
(in black lines). At ^ > 1000, the signal-to-noise does not vary 
substantially when considering different values of the param- 
eter c/, which represents the level of foreground subtraction. 



We will forecast parameter uncertainties in the following 
parameters tt = {R, crini?,6,r. Ay}. 

As a forecasting tool we will use a Fisher matrix anal- 
ysis, where the Fisher matrix is given by [47j : 



I 



dz 



{dC^^'\z)/d7r^)idCj''\z)/dTr,) 



(57) 



where /i and v run over the parameter modes. 

The rms uncertainty on the parameter tt^ is given by 
^{'^fi) — {^^jlY^'^ if the other parameters are marginal- 
ized. If the remaining parameters are assumed fixed, then 
the rms is a(7r^) = (F^^)"^/^ 

We consider a next generation polarization experiment 
with noise power spectrum given by Ap = 0.3/iK-arcmin 
and beam size Qfwhm = 1 arcmin. 

When assuming a LOFAR type experiment, the bias 
can be constrained at the 33% level, R at the 20% level, 
(j\nR at the 11% level. Ay at the 15% level, and r at the 
6% level, when the remained parameters are fixed. With 
a SKA type experiment, the bias can be constrained at 
the 11% level, R can be constrained at the 6% level, aini? 
at the 3% level. Ay at the 5% level, and r at the 2% level, 
when the remaining parameters are considered fixed. 

We show the error ellipses for the optical depth, the 
bias and the width of reionization in Figs. 11 and 12 



The Planck priors are shown in dashed lines. For a LO- 
FAR level of noise, and an experiment that measures the 
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FIG. 11: Forecasted uncertainties on the width of reionization 
parameter /\y and the optical depth r (assuming that the 
other parameters are fixed) for a LOFAR type experiment (in 
black lines) and a SKA type experiment (in red lines). For 
reference, the dashed lines correspond to the Planck priors on 
the optical depth. 
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FIG. 12: Forecasted uncertainties on the bubble bias b and 
the optical depth r (assuming that the other parameters are 
fixed) for a LOFAR type experiment (in black lines) and a 
SKA type experiment (in red lines). For reference, the dashed 
lines correspond to the Planck priors on the optical depth. 

full sky, the constraints on the optical depth are compa- 
rable to the expected ones from Planck. For a SKA type 
experiment, the constraints shrink by a factor of 3. 



VIII. DISCUSSION AND CONCLUSION 

We investigated the correlation between free electrons, 
traced by the optical depth r, and neutral hydrogen, 
traced through the emission of 21 cm photons, during the 
epoch of reionization. To compute the cross-correlation 
we used a simple model where patches of ionized gas are 
represented by spherical bubbles. The cross-correlation 



will depend on the presence or absence of these bubbles 
(the 1-bubble term) and the clustering of bubbles (the 2- 
bubble term). As expected, the cross-correlation is neg- 
ative on small scales, where it is dominated by the shot 
noise of the bubbles. On large scales, the 2-bubble term 
can render the correlation positive as long as the effec- 
tive bias 6eff is large. Small bubbles at a fixed neutral 
hydrogen fraction imply a small bubble bias, hence the 
2-bubble term has a suppressed (positive) amplitude. A 
larger correlation could be driven by the ionization frac- 
tion, but within a bubble merger scenario the smallest 
bubbles are expected at early times, when the ionization 
fraction is small and the total matter power spectrum is 
suppressed. 

The anti-correlation peak, set by the sum of the 1-and 
the 2-bubble terms, depends critically on the distribution 
of the ionized bubbles. Consequently, a measurement of 
the cross-correlation allows us to probe the parameters 
relevant for the reionization history. In principle a mea- 
surement of a positive correlation at early times, could 
set fairly stringent constraints on the bubble size and the 
bubble bias, providing a complementary probe into these 
EoR parameters. 

One major obstacle in measuring the 21 cm emission 
from the EoR are the large foregrounds at these frequen- 
cies. For the auto-correlation, any detection requires a 
careful removal of foregrounds, which typically results in 
the removal of the largest modes along the line of sight. 
The advantage of the cross-correlation is that foregrounds 
in the measurement of r are weakly correlated with those 
in the 21 cm field. Therefore, the cross-correlation is 
less sensitive to the detailed understanding of the fore- 
grounds. 

In this paper we have computed the signal to noise 
of the cross-correlation using the estimator for inhomo- 
geneous reionization r^rn proposed by Ref. [21 . In our 
computation there is very little contribution from any 
positive correlation at large scales coming from 2-bubble 
term, and most of the signal comes from the shot noise. 
Because a measurement of the optical depth gets most 
of its signal from the long wavelength mode along the 
line of sight, we left the 21 cm foregrounds as a noise 
term. Although the signal to noise per mode is small, 
the large number of modes allows for a detection when 
considering a next generation 21 cm experiment cross- 
correlated with a CMB experiment that measures the 
polarization 5-modes in most of the sky. We expect that 
around the time SKA observes a large part of the sky, 
CMB experiments will have improved to the level that 
we are able to reconstruct a map of r^rn- Although the 
auto-correlation of both maps will give significant insight 
into reionization, cross-correlating these maps will give a 
complementary probe. We find that a measurement of 
this cross-correlation with a detector noise level of a LO- 
FAR type experiment on the 21 cm side and noise level 
of a next generation polarization type experiment on the 
CMB side constrains the width of the ionization history 
at the 15% level, the bias parameter at the 33% level. 
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the bubble radius at the 20% level, and the width of the 
bubble distribution at the 11% level. The constraint on 
the optical depth (at the 6% level) is similar to what we 
expect to see with the Planck satellite. The previous re- 
sults are considering that both experiment measure the 
full sky. If they measure part of the sky, the constraints 

— 1/2 

will degrade as f^^^ . Finally, a SKA type experiment 
will be able to shrink these constraints by a factor of ^ 3. 
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Appendix A: Dependence on R and ainR 

The choice of a log-normal distribution is motivated 
in part by simulations by Refs. [32 and [29]. In this 
paper we have made the assumption that the distribu- 
tion is static (i.e. that cr\nR and the distribution it- 
self do not change in time). We have already seen that 
variations in (Ji^r change both the location and the am- 
plitude of the peak of the 21 cm-r cross-correlation. In 
this section we will consider implications of varying EoR 
parameters with a log-normal and a normal bubble dis- 
tribution on the l-and-2-bubble terms. In particular, we 
will compare both contributions to the amplitude of the 
anti-correlation peak. 



1. Log- normal distribution 

We first start by investigating the implications of a 
log-normal distribution for the bubble radius. 
The bubble distribution is given by 



1 1 



-[ln{R/R)f/{2af^^) 



(Al) 



R 



Given this distribution we can compute the average bub- 
ble size: 

(H) = J dRP{R)Vk{R) = l^e9<«/2 (A2) 

To address the dependence of the resulting correlation 
function, we also need the variance 



{Vi) = I dRP{R)Vi{R) 
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(A3) 



The amplitude of the 1- and 2-bubble terms, and the 
relevant scale where these peak, strongly depend on the 
window function 



Wnik) 



{kRY 



[sm{kR) - kR cos{kR)] ( A4) 



Recall the the volume averaged window function and 
window function squared (shown in Figs. 13 and 14) are 
defined as 

{WR){k) = r dRP{R)V,{R)WR{kR),{A5) 

\^b) Jo 



and 
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dRP{R)V^{R)Wl{kR) (A6) 



{V,)^ JO 

The total correlation function can be written as 
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^2b 



(A7) 



where the 1-bubble contribution consists of two relevant 
terms: the shot noise and the power spectrum P§§ given 
by Eq. (I36|. 



Let us define the following relevant ratios: 



TD^b—sn p26 / r>sn 



and 
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Ib—sn 



sn 
Xipy 



(A8) 



(A9) 



where is the contribution from the shot noise, which 



is given by 



-psn 
-^XiP 



-{x,-xl){Vt){Wl){k) 



(AlO) 



In the limit of small comoving momenta (large scales), 
(Wr) ^ 1, and {Wl) ^ {V,^)/{V,)\ 
For small scales, we have: 



{W^^){k,R,ainR) 



9 



(All) 



Roughly speaking, we know that the contribution from 
the shot noise term will be constant and have a peak at 
some characteristic scale after which it will decrease as 
(X 1/k^. Furthermore, up until that characteristic scale, 
the amplitude of the shot noise is boosted with respect 
to ah the other terms as {V^)/{Vi,)'^ = e^'^^^R. 

The characteristic scale is determined by equating the 
two limiting cases [48 , i.e. 



(A12) 



9 2k^ 
This tells us that the shot noise roughly peaks around 
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FIG. 13: {WrY for different values of 
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FIG. 14: {Wr) for different values of cr\r,R. Note that the 
maximum amplitude in the limit kR <^ 1 grows exponentially 
in cTini?, which will be relevant for the relative contribution of 
the shot noise with respect to the 1— and 2— bubble terms. 



The second contribution to the 1-bubble term comes 
from Pss{k)^ which is given by 



(fk' 
(2^ 



{Wi,){k')Pss{\k-k'\) (A14) 



As it was shown in Ref. [29 , in the large scale limit 
we have: 
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In this limit, {W^} {Vb)/{Vb)^, which allows us to put 
the following constraint on the amplitude of Pss{k) 
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FIG. 15: k^Px^p{k)/27v^ for = 1 Mpc, 6 = 6 at different 
values of ainR. 



The relative peak amplitude between the two contribu- 
tions in the 1-bubble term is therefore given by 



H = P66/Pxi; 



Jo 



k^dkPssik) (A17) 



Since kp^ak depends on the bubble radius and on aini?, 
so does the relative contribution. Generally speaking, a 
narrower distribution (with smaller Rq) leads to a larger 
contribution from Pss to the total 1-bubble term. That 
being said, even for very narrow distributions and very 
small average bubble radius we find that the total con- 
tribution to the peak does not exceed more then a few 
percent, i.e., in realistic scenarios the shot noise term 
dominates the total 1-bubble term (see e.g. Ref. [49 ). 
In Fig. 15 we plot the cross-correlation Xtp for different 
values ofoini^, confirming our estimate for the peak sale 
in Eq. (|A13[). 



There is one caveat, which is that Pss does not drop as 
fast as the shot noise, hence at small scales this term can 
contribute more. In fact, it is this term the responsible 
for the turnover at small scales of the r-r and 21-21 auto 



power spectra (see Figs. 16 and 17). 



We remind the reader that the 2-bubble contribution 
to the Xtp cross-correlation is given by: 



p26 



-xjj[\nxHb{WR)^lfPss{k) 
XH[\nxHb{WR)^l]Ps6{k) 



(A18) 



The minimum value P^^ is reached when xh = 
g- 1-1/6^ while for the shot noise term Xe = 0.5 repre- 
sents the peak value. We can now compare the 2-bubble 
terms to the shot noise: 

P5s{kpeak)e-^''^'-^'^ (A19) 



When considering R = 
that at z = 11, R'^b-sn 



1 Mpc and cti^r = 0.5 we find 
^0.6, which comes very close to 
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FIG. 17: k^P^^(k)l2i{^ for = 1 Mpc, 6 = 6 at different 
values of aini?. 



what we actually observe in Fig. jT] Note that the ratio 
depends on the redshift considered, but in general, we 
find that the 2-bubble term remains of the same order as 
the shot noise term. 

Can the 2-bubble term ever dominate over the 1-bubble 
term? Since the ratio goes as XjB? and decreases ex- 
ponentially in CTini?, for bubble distributions with small 
bubble radius and narrow width, we find that the 2- 
bubble term can easily dominate the total correlation 
function. In a more realistic scenario in which the ra- 
dius and the variance depend on redshift, we expect the 
2-bubble term to be increasingly important to the cor- 
relation function in the early stages of reionization, in 
other words, when the correlation function is dominated 
by points that live in two different bubbles. At the onset 
of reionization, the bubbles are small and it is more prob- 
able to find two points that live in two separate bubbles. 
As bubbles merge and grow, it becomes more likely that 
the correlation function has contribution from points that 
are in the same bubble. The authors in Ref. [6 make 
this distinction, and divide the reionization model in two 



regimes separated by the average ionization fraction. 

Concluding, we see that the location of the peak of 
the correlation function is roughly set by kp^ak^ given in 
Eq. (A13), and we note that the location of the peak is 
almost equivalent for the 1- and 2-bubble terms. 

We have shown that in the case of a log-normal dis- 
tribution the contribution from P§§ to the peak ampli- 
tude is small compared to the shot noise and as such, 
to the overall correlation. Since the shot noise grows 
as (X and exponentially in the width of the distri- 
bution <Jini?, we find that assuming smaller values for 
these parameters lead to rapidly increasing contribution 
of the 2-bubble term compared to the 1-bubble term. 
Since decreasing both of these parameters also increases 
the value of k^eak-, assuming a narrower distribution of 
bubbles with a smaller average radius results in a corre- 
lation function that is dominated by the 2-bubble term 
and peaks at smaller (physical) scales. These findings are 
consistent with the expectation that larger bubble imply 
a larger shot noise. 



2. Normal distribution 

Simulations show that bubbles are well traced by a 
log- normal distribution at early times 22| , while at later 
times the distribution can transition to a normal distri- 
bution. Since at late times, large radii dominate reioniza- 
tion, we expect the shot noise to become more dominant. 

A normal distribution is given by 



P{R) 



1 



(A20) 



As expected, for a Gaussian, all relevant quantities are 
much closer to the distribution values (e.g. R) deviating, 
by definition, at most 1 sigma. 

The average bubble volume is given by 



(A21) 



We will assume that a normal distribution is only valid 
for R > IMpc. In this limit, the above equality holds, 
even for a bound probability function, as long as < 1. 
This cutoff is consistent with observations. The expres- 
sion above can easily be understood by Wick expanding 
the 3-point function {R^). 

Similarly, for the variance we obtain: 



f (l5i?V| 



(A22) 



A gaussian distribution allows us to analytically com- 
pute the volume average window function: 



{WR){k) 



3e 



+ l) sin {kR) - kR cos (kR)] 



P {2,Ra\ 



(A23) 
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At large scales, 



The variance is given by: 



lim {Wr) = 1, 



while at small scales, 



lim (Wr) 



ScTpe 2^ sin 



3in (kR) 



kR (3cr| + l) 



(A24) 



(A25) 



ge-2k^al |g2/cV^ [^2 ^^2 ^ ^2^>) f^^P^j, + l) sin {2kR) + [e {R^ - 3al) - Ak^aj, - l] cos {2kR) } 



2k^ {3Raj^ + ' 



(A26) 



r 



At large scales, the variance becomes: 
lim (W^l) = (F,2)/(H>^ 

and at small scales: 



(A27) 



average does not give an accurate enough answer. There- 
fore, we will use the following best fit: 



(A31) 



lim (W^) = 9{R' + (ji,)/{2k\R' + 3^)^ (A28) 

By equating the two limiting cases, we can derive the 
peak scale for the variance of Wr: 



uvar 



2 (I5i?4cr| + 45i?2cr4 + i?6 + I5cr6 ) 



1/4 



(A29) 



The maximum peak value of the shot noise is then given 

by 



7,3 psn (h.var \ 



(A30) 



A similar approach for the peak scale of the window 



This and the derived maximum value of the shot noise 
immediately allow us to put a constraint on the ratio 
between the 2-bubble term and the shot noise, 



p26— sn p26 / psn 



Pss{k-/eak) 

'i-XH {v,){w^R){k;Zk) 

(A32) 



From this expression, we find that the shot noise is 
significantly larger than the 2-bubble term around the 
peak scale. This nicely fits into the previous picture, 
since the normal distribution of the bubbles is only phys- 
ical at late time in the reionization history when ^ > 1 
Mpc. At those times the 1-bubble term (dominated by 
the shot noise) should make up the largest contribution 
to the total cross-correlation spectrum. 
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